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Abstract 

We give an estimate of the quantum variance for d-regular graphs quantised with 
boundary scattering matrices that prohibit back-scattering. For families of graphs 
that are expanders, with few short cycles, our estimate leads to quantum ergodicity 
for these families of graphs. Our proof is based on a uniform control of an associated 
random walk on the bonds of the graph. We show that recent constructions of Ra¬ 
manujan graphs, and asymptotically almost surely, random d-regular graphs, satisfy 
the necessary conditions to conclude that quantum ergodicity holds. 


1 Introduction 

Quantum graphs have been suggested as an ideal model for problems in quantum chaos 
[1]. It is therefore surprising that there is no general theorem analogous to the Quantum 
Ergodicity theorem of Snirel’man, Zelditch and Colin de Verdiere [2, 3, 4] in the quantum 
graph setting. Quantum ergodicity has been proved for a special class of graphs derived 
from 1-dimensional maps [5] , and a general argument has been presented based on physical 
methods of supersymmetric field theory [6, 7]. On the other hand it has been proved that 
quantum ergodicity does not hold for quantum graphs with a star-like configuration [8], 
and entropy bounds—which control the extent to which eigenfunctions can localise—have 
been derived for a few different families of quantum graphs in [9]. A recent article proves 
quantum ergodicity for the somewhat related problem of eigenfunctions of the discrete 
Laplacian on combinatorial graphs [10]. 

Quantum ergodicity is one of the few universal results in quantum chaos. It implies a 
weakened form of the semi-classical eigenfunction hypothesis [11, 12], and can be stated 
in the following form: let (f n be an orthonormal basis of quantum wave-functions with 
energy levels E n , A an observable, with classical average A. Then 

liSo e } 22 !<*"'- A f= °- < L1 > 

provided that the classical dynamics are ergodic. (For the reader who prefers to keep a 
concrete example in mind, one can take <f> n to be a sequence of normalised eigenfunctions of 
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the Laplace-Beltrami operator on a compact Riemannian manifold of negative curvature, A 
can be a zeroth-order pseudo-differential operator with A the mean value of its principal 
symbol, and the classical dynamics are the geodesic flow on the manifold. Quantum 
ergodicity in this setting was proved in [2, 3, 4].) One implication of (1.1) is that one can 
extract a density-one subsequence 1 of wave-functions that becomes equidistributed in the 
large energy limit. This is equivalent to the semi-classical eigenfunction hypothesis for 
that subsequence [13]. 

Following the original proofs of quantum ergodicity, in the manifold setting, the result 
has been proved for a variety of situations, including Euclidian billiards [14, 15], quantised 
torus maps [16, 17, 18], and quantised Hamiltonian flows in M n [19]. 

In the present article we prove, deferring a precise statement of results to the following 
section, a quantum ergodicity theorem for quantum graphs quantised with the non-back- 
scattering boundary conditions introduced in [20], provided that the underlying graphs are 
expanders [21] and have not too many short cycles. These conditions are similar to those 
demanded for the proof of quantum ergodicity for eigenfunctions of the discrete Laplacian 
on combinatorial graphs in [10], although the method of proof there is quite different. 

2 Notation and statement of results 

In order to fix notations we briefly describe the main definitions that we will use in this 
work. For further background information on quantum graphs we refer the reader to the 
recent book [22]. 

2.1 Quantum graphs 

A quantum graph is a metric graph equipped with a differential operator acting in a space 
of functions defined on the bonds of the graph. We will denote by TJ and IB respectively 
the set of vertices and bonds of the graph, with |QJ| = n and |33| = B. The vector of 
bond lengths will be denoted L = (L b )b& b where each L b > 0. For us, all graphs will 
be undirected and simple, which means that no multiple edges are allowed, nor are loops 
connecting vertices to themselves. Furthermore, we will avoid considering bipartite graphs. 
This means that the vertex set cannot be partitioned into two sets with no connections 
within those sets. 

Our focus will be on d-regular graphs, d ^ 4, which are graphs where each vertex is 
connected to d other vertices. This imposes the (trivial) constraint 

2 B = nd. (2.1) 

Identifying each bond b of a graph Q with an interval [0, L b \ we can define spaces L 2 {G) 
as the direct product of interval L 2 -spaces. We will consider metric graphs acted on by 
the one-dimensional (positive) Laplace operator on intervals. The associated eigenvalue 
problem reads 

- = k 2 ^ b , b € ®, (2.2) 

and solutions are bond-wise waves, with vertices as scattering points. 

1 Meaning that the number of terms of the subsequence in a sufficiently large interval is asymptotic to 

the number of terms in the interval. 
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A vertex scattering matrix for a vertex v of degree d is a d x d unitary matrix a v . 
where the vector of amplitudes of incoming waves a in E C d is related to the vector a out of 
outgoing amplitudes by 

a out = a v & in . (2.3) 

By considering the 2 B directed bonds the bond-scattering matrix S' is a 2 B x 2 B matrix 
with 6cth entry equal to 0 if bond b does not feed into a vertex v that bond c leaves, and 
otherwise is the corresponding element of the matrix a v . The matrix S so-constructed 
is unitary, as a consequence of the unitarity of the cr v . The quantum evolution operator 
U = XJ{k ) is the 2 B x 2 B matrix whose fecth entry is 

U(k) bc = e ikL ”S bc , (2.4) 

where L b is the length 2 associated to the bond b. A common convention is to order the 
directed bonds so that bonds b = B + 1,..., 2 B are the reversals, in order, of the bonds 
b = 1,..., B. In that case, we can write U(k ) as the matrix product 

U(k) = D(k)S, where D(k) = l Q e ifcdiagL J , ( 2 - 5 ) 

however, we will sometimes adopt a different ordering for the directed bonds. 

We define the spectrum of the quantum graph to be the set of non-negative values k 
for which the condition 

det (U(k) - I 2B ) = 0 (2.6) 

is satisfied. We label the points in the spectrum as k m , m = 0, 1, 2,..., with multiplicity, 
ordered so that 

0 < k 0 < ki < k 2 < • • • (2.7) 

Condition (2.6) is equivalent to U(k m ) having an eigenvalue 1, and we define <P m E C 2B 
to be the corresponding eigenvector normalised so that ||^m||c 2B = or if = k m+ \ = 
■ ■ ■ is a multiple root of (2.6), take the corresponding ... to be an arbitrary 

orthonormal basis of the eigenspace at 1. 

If <7„ is a unitary matrix satisfying cr 2 = /^, then the procedure above is equivalent 
to choosing a certain self-adjoint extension of the Laplace operator, in the sense that 
the spectrum defined above is the eigenvalue set of the self-adjoint Laplace operator, 
with correct multiplicity (except possibly for the eigenvalue 0, see [23, section 5]), and 
the components of the vectors are the amplitudes for the wave solutions to (2.2), 
satisfying the boundary conditions implicitly specified by the choice of extension. To give 
an example, the boundary conditions: 

• V’b i (0) = i’bj (0) f° r every pair of bonds bj originating at the same vertex, and, 

• ^(0) = 0, where the sum is taken over all bonds originating at a vertex, 
leads, for a vertex of degree d, to the vertex scattering matrix cr v , where 

()ij — ~i dij . 

2 Directed bonds have the same length as their undirected counterparts. 


3 



An alternative point of view treats vertices as scattering centres, where unitarity is 
a necessary condition on the matrix of transition amplitudes to ensure probability con¬ 
servation. This point of view legitimises the assignment of arbitrary unitary matrices to 
vertices [24, 25], which allows for greater flexibility by choosing matrices with advanta¬ 
geous properties. In [20] a class of scattering matrices were introduced with the properties 
that all diagonal entries are 0, and all off-diagonal entries have equal complex amplitude. 
These were referred to as equi-transmitting matrices. 

The diagonal elements of a scattering matrix give the reflection amplitude for a wave 
to be back-scattered into the reversal of the original bond. By setting these elements to 
zero, back-scattering is prohibited in the corresponding quantum graphs. In [20] it was 
speculated that these equi-transmitting quantum graphs would lead to new advances in 
the study of quantum chaos on graphs. Our proof in the present article of a quantum 
ergodicity theorem for quantum graphs with equi-transmitting boundary conditions can 
be considered as such an example. 

Equi-transmitting matrices of size d x d have been proved to exist in [20] for d = 2 n , 
d = P+1 where P is an odd prime number, and any d for which a skew-Hadamard matrix 3 
exists. Equi-transmitting matrices are further known to exist for d = P n + 1 for any n > 1 
and P any prime [26]. 

2.2 Notions of graph theory 

We will appeal to a few notions of graph theory, collected here for reference. Classical 
graph theory is concerned with combinatorial graphs, i.e. without reference to any bond 
lengths. The connections are encoded by a n x n matrix C called the connectivity matrix , 
whose ijth entry is 1 if vertices i and j are connected, and 0 otherwise. If the graph is 
without multiple edges or loops, and d-regular, then C is symmetric, and each row contains 
precisely d l’s. 

The (combinatorial) spectrum of a graph can be defined in a few ways, which are 
coincident if the graph is regular. We shall define it as the set of eigenvalues /xi, of 

C. As C is symmetric, the eigenvalues are real, and we order them in decreasing order so 
that 

- d < n n ^ < hi = d. (2.8) 

The multiplicity of the eigenvalue d is the number of connected components of the graph 
(so that /ii = d —every graph has at least one component), and fj, n = —d if and only if 
the graph is bipartite. The eigenvalues of C excluding +d will be called the non-trivial 
spectrum. 

We will consider sequences of graphs indexed by an increasing number of vertices 
n —>• oo. Such a sequence of d-regular graphs is called a family of expanders [21] if there 
exists a constant /3 > 0 such that the non-trivial spectrum of each graph in the sequence 
is contained in the interval [—d + (3, d — (3]. If we can take (3 = d — 2 y/d — 1 the graphs are 
called Ramanujan. Ramanujan graphs are extremal in this sense, since for an increasing 
sequence of graphs the Alon-Boppana bound [27, 21, Theorem 5.3] implies 

liminf /i 2 ^ 2 \/d — 1. (2.9) 

71—> OO 

We shall refer to cycles on a graph, which are closed paths without back-tracking. We 
define the set Cs,t to be the set of bonds b E 21 of a graph that lie on a cycle of length at 

3 A skew-Hadamard matrix H is a matrix whose entries are ±1, with orthogonal columns, and satisfying 
H + H T = 2I d . 
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Figure 1: An illustrative example for the definition of 7b, t- The bond bo belongs to 7b, t 
for each t ^ 4. 

most t. The girth of a graph is the length of the shortest cycle. In particular, this means 
that CB,t = 0 whenever t is less than the girth. 

We also define the set 7b, t to be the set of directed bonds bo such that there exists 
t\, t 2 with ti +12 = t and bond bo is a distance at most t,\ from a cycle of length at most 
2t 2 (see figure 1). 

The sets 7b, t and Cs,t both give a measure of the number of short cycles. The set 7b, t 
is more useful for our purposes, but C B ,t is easier to understand. Fortunately, their sizes 
are related, as the following lemma (which does not give the sharpest possible statement) 
makes clear. 

Lemma 2.1. Consider the sets Cs,t and 7 b, t defined above, for a d-regular graph. Then 

(d — lT 1 

ITbA < 1 d _’ 2 \C B ,2tl ( 2 - 10 ) 

where \ ■ \ denotes the number of elements of a set. 

Proof. Clearly, 

\r B ,t\ < ^2 (d - l) tl |Cs,2t 2 |- (2-11) 

tl+t 2 =t 

Furthermore, t 2 2 for graphs without loops or multiple edges, so we have 

t 

\r B , t \< E( d -i) t_t 2 l c ^ 2 | 

i 2 =2 

t 

\C B ,2t\ ^(d- l)*”* 2 , (2.12) 

t2= 2 

since \C B , 2 t 2 \ ^ \C B , 2 t\- We sum the geometric series to get 

I?^^^[((d-i)*" 1 - 1 ) 

(d — l) t_1 . 

< d _ 2 \C B ,2t 1- (2-13) 

□ 
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2.3 Quantum Ergodicity for quantum graphs 

Observables on a quantum graph will be functions that are constant on directed bonds, 
which can be represented by members of C 2B . For such an observable / € C 2B , the 
quantisation of /, denoted Op (/) is simply the diagonal 2 B x 2 B matrix containing the 
entries of /: 

Op(/) :=diag{/}. (2.14) 

Let {(f>j(k)} 2B 1 be an orthonormal basis of eigenvectors of the matrix U = U(k). We 
define the quantum variance as, 

\ \ r K 2B 

V ^ ,B ^ := 2BK m ooK / 2 |( < M fc )’°P(/)^#))c 2S - 2^ T rOp(/)|^ dk. (2.15) 

wo ° J o j=1 

While V(f,B) defined by (2.15) does not seem immediately analogous to (1.1), it was 
proved in [28] that averaging the second moment appearing in (2.15), involving the eigen¬ 
vectors of the matrix U(k), over a large window of k values is equivalent to averaging an 
expression similar to (1.1) involving the eigenstates <Fq, ^i, $ 2 , ■ ■ at least for graphs with 
incommensurate bond lengths. We take (2.15) as the starting point of our investigation. 

It turns out (see [5] for example) that one cannot expect that V(f,B ) = 0 for any 
individual fixed graph. For a fixed graph with Kirchhoff boundary conditions, a complete 
classification of limiting measures induced by subsequences of eigenfunctions has recently 
been obtained in [29]. Consequently, we consider a family of graphs, indexed by a sequence 
of increasing number of bonds B —>• 00 . Quantum ergodicity (sometimes called asymptotic 
quantum ergodicity) for quantum graphs means that H(/, F>) —^ 0 as .B —>■ 00 along this 
sequence, [5, 7]. 

Our main result, stated below, estimates the size of the quantum variance for graphs 
quantised with equi-transmitting matrices, in terms of certain graph-theoretic properties. 


Theorem 2.2. For d > 3 consider a d-regular connected, non-bipartite, simple graph, 
on B bonds, quantised with equi-transmitting scattering matrices. Let f € C 2B be an 
observable, satisfying |/&| < n for all b, for some k > 0. Let T > 0, and suppose that all 
non-trivial eigenvalues /j,; of the connectivity matrix of the graph satisfy \ni\ ^ d — fd for 
some (3 > 0. Then the quantum variance (2.15) satisfies the main estimate, 


V(f,B ) 



K 2 (d — 1) T |Cb,2t| 
BT 2 


(2.16) 


where Cb, 2 T denotes the number of cycles as was described in section 2.2 above. 


It is natural to impose a uniform boundedness condition on observables for which we 
wish to prove that quantum ergodicity to hold. This means that the parameter k in the 
statement of theorem 2.2 is an absolute constant, independent of B. However, the estimate 
(2.16) makes it clear how to consider observables that grow mildly as B —» 00 . 

Theorem 2.2 proves quantum ergodicity for families of graphs for which a parameter 
T —>• 00 can be found as B —)■ 00 , in such a way that the quantities on the right-hand 
side of (2.16) become negligible. We give two examples for which this is the case: the 
families of Ramanujan graphs constructed in [30], and random d-regular graphs [31, 32]. 
We discuss these examples further in section 6 below. 

We also mention that theorem 2.2 does not require that the bond lengths of the graph 
are linearly independent over Q, in contrast to many other results in this field, although 
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if the bond lengths are incommensurate, some values of constants can be improved—see 
the comment at the end of subsection 6.2. 

We remark that we excluded bipartite graphs from our consideration in the introduc¬ 
tion. Our methods could be extended to include bipartite graphs, but there the notions of 
ergodicity would need to be generalised, since bipartite graphs can support non-uniform 
invariant states. 

All results in the present paper would hold in the case d = 3, were it not for the fact 
that no equi-transmitting matrices of size 3x3 can exist (as a short calculation shows). 
Therefore theorem 2.2 is stated with condition d > 3, although in later parts of this work 
d can be set equal to 3. 

In most quantum ergodicity results, such as [2, 3, 4], the proof naturally separates into 
a semi-classical part, and a dynamical part. In the semi-classical part, a correspondence 
is established between quantum and classical time evolution, and in the dynamical part, 
ergodic properties of the classical dynamics are invoked to prove quantum ergodicity. We 
shall present our proof of theorem 2.2 in this manner. In section 3 below, we relate (see 
proposition 3.3) the quantum variance to a classical random walk on the bonds of the 
graph. In section 4 below, we analyse the equidistribution of the random walk, proving a 
uniform decay estimate that allows us to obtain (2.16). The final steps of the proof are 
carried out in section 5. 


3 Semi-classical argument 


Our aim is to prove, for a suitable class of /, that V(f,B)->0asB—>- oo. Without loss 
of generality we can and will assume that / is chosen so that Tr Op(/) = 0, which will 
simplify the notation somewhat. 

Our main estimate for the quantum variance goes back to an idea from [33] (see also 
[18] for a similar application of this idea). Let T > 0 and 


w T (t) := 



1*1 < T, 
otherwise. 


(3.1) 


Then for any 2 B x 2 B unitary matrix U, and any orthonormal basis of eigenvectors 

of U, 

1 2B t 

— ^|(^,Op(/)^) c2B | 2 £ ^t(*) Tr(Op(/)*t/*Op(/)t/ _t ). (3.2) 

3=1 t=—T 

(In order to make this paper as self-contained as possible, we include a proof of (3.2) in 
an appendix. See lemma A.l). Therefore, we can estimate the quantum variance by 

E ^t(*)( Jim^- / Tr(Op(f)*U(k) t Op(f)U(k)~ t ) dkj . (3.3) 

n t=—T _> °° "'° 


The limit in (3.3) exists as the integrand is an almost-periodic function of k. 

Let us define the matrix M® for t € No as the 2 B x 2 B matrix whose be th entry is 

|c,w - |2dt - (3 - 4) 

i.e. the average of the square of the be th element of the matrix U ( k) 1 . Then we can rewrite 
the quantity inside the brackets in (3.3): 
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Lemma 3.1. Let Op(/) = diag{/} where f G C 2B and let M® be defined as above. Then 

lim 1 [ K Tr(Op(f)*U(k) t Op(f)U(k)~ t ) dk = (/.P/)^, (3.5) 

a-kx. A Jo 

/or tfi 0. 

Proof. We denote A = Op(/). Expanding the trace, we have 

Tr(F*U t FU~ t ) = ^2 F bobo U bobl U blb2 • • • • • • ^ l6o , (3.6) 

b 0 ,—,b t -i 

Co,...,C(_l 

where the multi-sum runs over all possible choices of 2B bonds for each of bo, ..., bt-i and 
Co,..., ct- 1 - We have also the expansions 

U bo c 0 — 5L U b0bl U bl b 2 ---U bt _ic 0 i (3-7) 

bi,...,bt~i 

and 

Uboco = U b0Ct _iUct-,ct-, ' ' ' U C1C0 

= E (3.8) 

Cl,...,Ct-l 

so that 

J K MF*U(k) t FU(k)- t )dk= ^ A 6ofto (^hml J K \U{k)\ QC0 \ 2 dfc) A CoCo 

° f>0,CO 0 

= </,M (4) /) c 2B , (3.9) 

□ 

The matrix is not easy to work with, so we introduce a second matrix M where 

M bc ■= \U bc \ 2 . (3.10) 

Because U is a unitary matrix, the matrices M and M® are both doubly stochastic. If 
the graph does not contain too many cycles, then the matrix M® is close to M t , in the 
following sense: 

Proposition 3.2. Let T G N and suppose that f G C 2B satisfies the bound \f b \ ^ n, 
6 = 1,..., 2 B for some k > 0. Provided that the graph is quantised with scattering matrices 
that prohibit back-scattering, then 

2 k 2 

(/, M (f) /) C 2 s - (/, M*/)c 2S ^ (d-2)(d- !) ~~ 1 )*I C b.2t|, (3.11) 


/or all t = 1, 2,..., T. 




Figure 2: An example of a path and its return with t = 4. In this example the return path 
is different to the outward path. If the graph does not contain too many cycles of length 
at most 6 then this happens only rarely. 




Figure 3: Two different possibilities for excursions along a path from bo to Co- In (a) there 
is back-scattering; in (b) no back-scattering occurs but we note that the bond bo belongs 
to the set Tb, 5 - 


Proof. By (2.5), we can expand \U(k)^ \ 2 as 


I U{k) 


t |2 
b 0 c 0 I 


S bobl S blb2 ■ ■ ■ S bt _ 1CQ e lk<<Lbl +Lb2+ '" +Lb t~i +Lc o ) 


SC]_CQ SC2C1 


' S boCt _ 1 e 


— ifc(L C Q+L cl H-) 


(3.12) 


Because S bc = 0 if directed bonds b and c are not connected, we can think of the right-hand 
side of equation (3.12) as a weighted sum over paths connecting bo to Co and a return path 
(see figure 2). 

The crucial step in our argument is to demonstrate that, with few exceptions, the 
return path goes back over the same bonds in the reverse order as the outward path. This 
might not happen if, along the path, there are places where at least two distinct excursions 
from the same vertex are made, in the sense that removing the excursions gives a shorter 
path from bo to Co (see figure 3). Since the graphs are quantised without back-scattering, 
the excursions may not consist of self-retracing sections, so the only possibility is if the 
excursions contain short cycles. 

A second mechanism whereby the return path may differ from the outward path is if 
the path contains a short cycle of an even number of steps where the outward path takes 
one route, and the return path takes a different route, as could happen in the situation 
depicted in figure 2. 

If the bond bo does not belong to the set Tb,*, then there are no cycles close enough 
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to bo to allow either possibility. In this case, the return path has to be the reversal of the 
outward path. This means that c\ = bt- 1, C2 = h-2, ■ • •, <h- 1 = b\ and so 


Mk)l oc /= J2 I 5 w>iI 2 |Sm 2 I 2 ---I‘V 1 c 0 | 2 

= (MXco, (3.13) 

since 1S& c | = |[/& c |. Since (3.13) is independent of k, the average in k in (3.4) has no effect, 
and we have 

(mW) =(M*)w (3.14) 

\ / O 0 C 0 

for bo ^ 7 b, t- Therefore, let us define a matrix Rw by 

R (t) = jQ-(t) _ m*, (3.15) 

and let us consider (R^ f)b, the 6th component of R^f. We have proved that 

{R (t) f) b = 0, (3.16) 


if b ^ 7 b, t- If b E 7 b, t then we can be sure that 

ICR (t) /)&l<2«, (3.17) 

since the matrices M^) and M (and hence M f ') are doubly stochastic. Hence 






^ 2k 2 |Tb^| 


(3.18) 


using lemma 2.1, 


‘Ik 2 

^^2 ( d- 1)t - 1 |C B) 2T|, (3.19) 

using the fact that \CB, 2 t\ \Cb,2t\ for t ^ T. □ 

As we shall see below, quantum ergodicity essentially follows if we can prove that 
(/, M^/) C 2 s = o(l) as t —>• oo. However, we are able to prove this decay for the matrix 
M t only, and proposition 3.2 provides the requisite link. We remark that this procedure 
is reminiscent of the recent proof of quantum ergodicity for ray-splitting billiards [15]. 
In that work the authors consider a probabilistic random walk on families of trajectories 
with the same end-points, with transition weights given to individual trajectories, and 
summed over families of trajectories; the same distinction as between our matrices M t 
and M^\ with quantum ergodicity likewise following from ergodicity of the latter class 
of random walk. In [15] the multiple trajectories arise as a result of splitting trajectories 
at the interface between one-or-more different billiard media; in our work the multiple 
“trajectories” arise due to the connectivity of the graphs. 

The main result of this section is as follows. 
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Proposition 3.3. Consider a d-regular graph on B bonds, quantised without back-scat¬ 
tering. For any T > 0, and any f € C 2B satisfying \ff\ ^ k for each b and Op(/) = 
diag{/} satisfying Tr Op(/) = 0, 

T 

V(f, B) ^ ^ Tr Op(/) 2 + \ fl ™T{t) (</, MV) C 2 b + O ( n 2 (d - 1 Y\C b ,2t\) ) • (3.20) 

n t =l 

Proof. Because of cyclic invariance of trace, and the symmetry of wt, we can write (3.3) 
(extracting the t = 0 term) as, 

V(f, B) ^ Tr °P ( /)2 + -^ ^(t) ( lim 1 f Tr(Op(/)*C/(A;) t Op(/)C/(fe)- t ) d/c) . 

-o t=1 wo ° io 

(3- 21) 

We then use lemma 3.1 and proposition 3.2 to get (3.20). □ 

To prove our quantum ergodicity result, we will need to understand the behaviour of 
M t f , which represents a random walk on the bonds of the graph. This we do in the next 
section. 


4 Dynamical argument 


4.1 Classical dynamics on a quantum graph 

The classical analogue of the quantum evolution is the Markov process on the directed 
bonds of the graph with the transfer matrix M, which is doubly-stochastic [34]. 

Since M is not necessarily normal, it will be convenient to work, rather than with 
eigenvectors of M, with its singular vectors, defined to be the eigenvectors of the symmetric 
matrix M T M. If we order the directed bonds in groups of d bonds departing from each 
vertex, the matrix M T M decomposes into block-diagonal form 


J = 


(d-iy 



( J 

0 ••• 

0 ^ 


— 

0 

J ••• 

0 

1 


V o 

0 ••• 

J 1 


d matrix 




/ d 

- i 

d- 2 


d -2 \ 

d 

- 2 

d- 1 


d- 2 

V d 

- 2 

d- 2 


d — 1 ) 


The eigenvalues of J are 1 with multiplicity one, and (d— 1)~ 2 with multiplicity d— 1, and 
the simple eigenspace is spanned by the vector (1,..., 1) T . Therefore the singular values 
of the matrix M are: 1 with multiplicity n, and (d — l) -1 with multiplicity n(d — 1). A 
basis for the eigenspace of M T M with eigenvalue 1 is given by the set {e±, ..., e n } C C 2B , 
where the j th component of e v is defined to be 1 if directed bond j points outwards from 
vertex v , and 0 otherwise (see figure 4). 

We could equally-well consider MM T , which has identical spectrum to M T M, and 
a basis of eigenvectors corresponding to the singular value 1 can be chosen to have zero 
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Figure 4: Support of the singular vectors e 3 and e.j associated to a vertex v. See main 
text for the definitions. 


components except for the incoming bonds of vertex v where the component is 1. We will 
denote these vectors by e v € C 2B for v = 1,..., n (see figure 4). 

Observables that are linear combinations of e±,... ,e n will be called evenly distributed 
around vertices in the following subsection. 


Proposition 4.1. If the scattering matrices on the quantum graph are equi-transmitting, 
the action of M on the vectors e v and e v is as follows: 


Me v = e 


VI 


Me,, = 


1 

d- 1 



(4.1) 

(4.2) 


where is a sum over vertices w connected to v. 


Proof. The action of M on a vector supported on a single directed bond b is to allocate 
an equal fraction to the d — 1 directed bonds feeding into the origin of b (not including 
the reversal of b ). For the vector Me v , each bond directed towards v gets d— 1 times the 
fraction l/(d — 1) of the weight 1 on each outward pointing bond in e v . The result is a 
vector of weight 1 on each inward pointing vertex to v. Hence Me v = e v . 

In a similar way, it is clear that Me v has weight 1 /(d — 1) on each bond directed 
towards a neighbour of v, except the outward pointing bonds from v (see figure 5). Such 
a vector can be written 



□ 


4.2 Observables evenly distributed around vertices 

Let 

Gi := spanjei, ... ,e n j C C 2B , 


(4.3) 
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Figure 5: The support of Me v . The vertex v is at the centre of the illustration. 


and let p : G\ —>• C n be the natural isomorphism. Let us also define (p \ G\C 2n by, 

f ai \ 


(p \ a>\e.\ + • • • + a n e n K y 


O-n 

0 


V 0 J 

We shall also work with the larger space G \, where 

Gi := span{ei,.., e n , e 1} ... ,e n } C C 2B , 
and define ip : C 2n —>• G\ by 


where 


ip ■ ( b ) aiei H-b a n e n + b\e\ H-b b n b n , 


a = 


a i 


V 


We note that ip fails to be invertible since 


ip I ~ 1=0, where e = 

—e 


/ 1 
V 1 


Lemma 4.2. If x = ( " ) € C 2n then 
' Y 


Proof. Let us define 


x)|| C 2s ^ \/d(||x|| c « + ||y||c")• 


r l = J2 x i e j 

3 = 1 
n 

3 =1 


(4.4) 


(4.5) 

(4.6) 

(4.7) 


(4.8) 


(4.9) 
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so that 



IIV’(®)llc 2B — Wv + vWcib — M ic 2S + 7 ? c 2S + 2 fKc ( 77 , fj) C 2 B. 

(4.10) 

Due to orthogonality of (ei,..., e n }, 



n n 

iMic 2 -® = y " i*j-i imcm = i y\ x j\ = ^ii x iic n ! 

3 = 1 i =1 

(4.11) 

and similarly 

C 2S = % €"• 

(4.12) 

We also have 

n 



(V,V) C 2B = Y, X iVj( e i^j)C 2B - 

(4.13) 


i,j = 1 

Now, since ej is supported on outward pointing bonds from vertex j and e t is supported 
on inward pointing bonds to vertex i, the only way that ( ei,ej)<£ 2 B can be non-zero is if i 
connects to j. We have, in fact, 


so 


(ei,ej) C 2B 


1, if * ~ J, \ = c 
0, otherwise, J 


(??, fj) C 2 B = ^ XiVjCij 

i,j=l 

= ( x , C y)c n - 


(4.14) 


(4.15) 


We get, 


\\-ip(x)\\ 2 C 2 B = d\\x\\ln + d||y||c" + 2 fRc(x,C'y) t > 
^ ^ll x llc n + d||y||c" + 2||x||c« HCyllc" 
^ ^ll x llc n + ^lly|lc" + 2d||x||c" ||y||c™ 
= d(||x||o + Hylic") 2 , 


proving (4.9). 

Let C be the 2 n x 2 n matrix, 


C 


0 d- i A 
^ £ T C ) ' 


Proposition 4.3. Let f E G\ and t = 0,1, 2,... Then 

i> (&<ptf)) = M f f. 
Proof. Let x € C 2ra with g = We first prove that 


= Mg. 

Indeed, suppose that 


g = x\ei H-h x n e n + xiei H-h x n e n , 


(4.16) 

□ 

(4.17) 


(4.18) 


(4.19) 


(4.20) 
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so that 


Mg = xiei 4- x n e n + L- ( ^ ej — e^j -\ -f 4 “ e?l )’ ( 4 - 21 ) 

jr^J 1 jf~n 

by proposition 4.1. However, with x = (xi,... ,x n ) T , x = (xi,... ,x n ) T and x = (x,x) T , 
we have 

c^( x ;4). (4.22) 


from which we see that if>(Cx) = Mg. 

To prove (4.18) in the case t = 0, it is immediate to observe that if{(p{f)) = /, from 
the definitions of if and (p. If we assume that ip(C t ip(f)) = M t f, then applying (4.19) 
with x = C t tp(f) proves (4.18) for the t + 1 case. □ 

Therefore, to understand the action of M on /, we need to understand the iterates of 
C. 


Let x be an eigenvector of C with eigenvalue g. Then 



Denoting 


So 



Zt.X 


z t -1 

d-l x 


(yt -1 + 


v-zt-l 
d -1 



yt 


Zt -1 

d- 1’ 


(4.23) 


(4.24) 

(4.25) 


(4.26) 


and zt = zt(g) is the solution to the recurrence 


_ gz t -i _ z t - 2 
Zt d — 1 d — 1 


(4.27) 


with initial conditions zq = 0 and Z\ = 1 . 


Proposition 4.4. If \g\ ^ d — /3 < d, then the solutions to (4.27) satisfy 


\zt{g) | 


< t 


( d- 1-/3 
V d — 1 


t-i 


(4.28) 


fort = 1 , 2 ,... 

Proof. We first consider the case |/x| 4 2 \/d — 1. Standard methods yield the solution to 
(4.27) in this case to be given by 


\Jd — 1 ^ u + y/u 2 — l\ ^ cj — \/w 2 — 1^ ^ 

2\iuf 2 — l y y Vd — f y y \Jd — l y y 


(4.29) 


15 



















where u> = \n{d — 1) 1 / 2 / ±1. Since 


y yfd^l 

we may apply the inequality 


u + y/u 2 — 1 \ ( u> — y/eo 2 — 1 \ y/u 2 — 1 


Vd — 1 I y/d — 1 


^-6* 


a — b 


^ t max{|a|, 16|}* 1 , 


to (4.29), to get 


|zt| ^ t max 




+ \/w 2 — 1 


y/d^l 


to 


— y/to 2 — 1 


For 1 < w ^ i(d — /3)(d — l) 1 / 2 we have 


w 


± \/w 2 - 1 


yfd ^I 




CJ 


yfd^l 


+ \/w 2 — 1 


t-1 


(d-l) 1 / 2 ' 


For such values of to, 


0 < uA - 1 < 


(d - P - 2) 2 - 4/3 (d-p-2) 2 


< 


4(d — 1) 4(d — 1) 


so that 


(4.30) 


(4.31) 


(4.32) 


(4.33) 


(4.34) 


u; 


db Vto 2 - 1 


yfd^l 


„ Id — /3 d- 2-/3 
< 2d- 1 + 2(d — 1) 

d — 1 — /3 


d- 1 ' 

A similar argument holds if — ^(d — /3)(d — l) 1 / 2 ^ to < —1. 
If M < 1 then 

w ± \/w 2 — 1 


1 ^ d - 1-/3 

\/d^I ^ d — 1 


Finally, in the case |/x| = 2\/d — 1, directly solving (4.27), we find 

t 


\zt(f*)\ = 


{d - l)(*-i)/ 2 ' 


Our main result of this subsection is the following: 


(4.35) 


(4.36) 

(4.37) 
□ 


Proposition 4.5. Let f £ G± with TrOp(/) = 0. If all non-trivial eigenvalues [i of C 
satisfy the bound \fi\ < d — (3 then we have 


\\M t f\\ C 2B ^ 2||/|| C 2Bt 


d- 1-/3 
d — 1 


t -1 


(4.38) 


fort = 1,2,3,... 
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Proof. Let x i,...,x n be an orthonormal basis of eigenvectors of C, where x\ = n l / 2 e. 
Since TrOp(/) = 0 we have (</>(/), e)c« = 0 and we can write 


<p{f) = a 2*2 H-h a n x n , 

where 

M/)llc" = l a 2h + • • • + \&n\ 2 ■ 

It is also easy to see that 

Wf\\&B=dMf)\\&. 

We have 

vif) = a 2 o") -1 - l " Qn | o) G ^ 2n • 

So, if i = 1,2,3,..., 

At _ 012 ( ~ Zt ~ i(^2)x 2 \ . «n / ~ Z t-lipnVn \ 

° d _ ! ^ (d - l)z t Qx 2 )x 2 l)^(A^)xn 7 ’ 

using (4.24). We now use lemma 4.2 to get 


\\M t f\\ C 2B= V’ C'V 


C 2B 


< 


€ 




a. 


1 2 I Zt.—l (/b ) 1^ 

' (d-1) 2 


1/2 


1/2' 


+ 5^M 2 M/b 

\j=2 


ME 


\j =2 

\ 1/2 

n 2 I (i-l)(d-l-/9) t_2 + i(d-l-)8) t_1 


Qi 


■ 1=2 


(d-1) 


t-i 


using proposition 4.4, 

^ 2||/|| C 2flt 


d- 1-/3 
d- 1 


t-i 


(4.39) 

(4.40) 

(4.41) 

(4.42) 

(4.43) 


(4.44) 

□ 


4.3 General mean-zero observables 


We would like to consider a more general class of observables than those belonging to the 
spaces G\. Let G 2 = Gi . Thus, G 2 is the space of singular vectors of M with singular 
value -py, and it follows that M acts on G 2 by contraction: 

Lemma 4.6. Let g £ Gj. Then, 

\\Mg\\cu> = 


Proof. Take g £ G 2 . Then, we have 

WMgW^B = (Mg,Mg) C 2B 


(g,M T Mg) C 2B 

Ml hs 

(d-if 


Finally, let us put the results of proposition 4.5 and lemma 4.6 together. 


(4.45) 

(4.46) 

□ 
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Theorem 4.7. Let f £ C 2B with TrOp(/) = 0. Assume that all non-trivial eigenvalues p 
of the connectivity matrix satisfy \p\ ^ d — /3 for some j3 > 0. Then there exists a constant 
K > 0 (that does not depend on t or B) such that for t^l, 

l|M7llc?» < A-||/|| c ,»t (4.47) 

Proof. The conditions of the theorem guarantee that 

l|M+|| c „ « ( ‘ J 1 (4.48) 

for any g € G\ for j € N, due to proposition 4.5. 

We decompose / according to C 2B = G\ © G 2 as / = /o,i + fo, 2 , where /o,i £ G 1 and 
/o ,2 £ G 2 , and inductively defining the sequences {./"y. 1 }y^=o an( 4 {/j,2}j*Lo by 

Mfj, 2 = /j+ 1,1 + (4.49) 


where fj^ € Gi for j ^ 0, i £ {1, 2}. So upon each iteration, the component of fjp that 
does not remain in G 2 becomes fj+ 1,1 € G\. 

We have, by lemma 4.6 


WfoWh* 

(d-l)2 


I|m/,- 2 ||2 2S 


H/7+1,1 


112 

llc 2S 


+ 11/7+1,2 


112 

ll<C 2B 


5 s II/ 7 +I,* 


II 2 

IIC 2S ’ 


for i £ {1,2}, so, inductively, 


(4.50) 


Tj, 2\\C 2B 




I l/o,21 lc 2S < ll/llc 2S 


(d-iy (d-iy 


and 


IU 11 / II/7-1,2||c 2S ^ ll/llc 2S 


Acting on / iteratively, we have 


M f f = M t f o,i + MVo,2 

= MVo,l + M 4 " 1 /!,! + M*-7l,2 


t-1 


Af* ifj t 1 + f t , 1 + f t) 2- 
j =0 


(4.51) 


(4.52) 


(4.53) 
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Thus, we have, using (4.51) and (4.52), and (4.48), 


t—i 


\\M f f\\ C 2B ^ ^2 ll M< j fj, lllc 2S + \\ft,l llc 2B + II/ 1,2 I 
3 =0 

/I t-j-1 


^2 


3=0 
t -1 


d-1-0 
d- 1 


ll/j,i|lc 2fl 


2 II/IIc 2b 

(d-l)* 


(d ~ 1 — j ll/llc 2 - 0 . 2 ll/llc 2S 
* ^ d-1 1 (d-iy + (d-i)< 


3=0 


2 II/II C »_ r (d - 1 - py-' + mar = t-j, 


(d- 1) 


t-1 


r=l 


(d-ir 


^ 2t||/|| C 2B (d — 1) 

C (d- if (d-2-P) (d-lf 

€ 21 HZ|| C 2a(d — 1) ^ d-l-^ y | 2||/|| C 2 b 


(d — 2 — /3) V d-1 / ' (d-l)* 
Finally (to combine the two terms into a single), noting that 

t(d — l)(d — 1 — 13)* 


d-2-/3 


>4, 


for 1 ^ 1, we get 


for 1 > 1. 




2(d - 2 - /3) 


d — 1 


(4.54) 

(4.55) 

(4.56) 

□ 


5 Quantum ergodicity for equitransmitting expander graphs 

We are now able to prove theorem 2.2. To begin with, we need a few results on certain 
summations. 

Lemma 5.1. Let 9^1 and T > 0. Then 


£ w ‘ 

t= l 


T9 t+2 + 0- (T+1)9 t+1 

(^i? : 


(5.1) 


and consequently, if \9\ < 1, 


£»■ 

t=i 


0 

(^F 


(5- 2 ) 


Proof. The sum appearing in (5.1) is of a standard type (see formula 0.113 of [35], or 
[36, p. 33] for a derivation). Equation (5.2) follows by letting T —>• oo. □ 


Lemma 5.2. Let T > 0 and 


w T (t) := 



|f| < T, 
otherwise, 


(5.3) 
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and 9^1. Then 


J2 9t ™ T (t) = 


t =i 


2^2 


T - 1 + 0 T - T9 


Proof. We reverse the order of summation, to get 
T T—l 

Y O^xit) = Y 0 T ~ k w T (T - k), 


via k = T — t, 


t =l 


k =o 


e' 


T T ~ 1 




—k 


k =0 
T—l 


1 - 


T — k 


Y he 


-fc 


fc =0 


To evaluate the sum in (5.5) we use lemma 5.1. The result is 


t=i 


( 0- 1 - l ) 2 

9 {T — l + 9 T — T9' 


J 1 2 


(1-0) 2 


(5.4) 


(5.5) 


(5.6) 


□ 


Proof of theorem 2.2. We recall equation (3.20) which provides the main estimate for 


1 1 1 

V(f, B) ^ — Tr Op(/) 2 + - Y Mt) (</> M t /) C 2 B + O ( K 2 (d - 1)*|C B , 2T |) ). (5.7) 


t=i 

In order to estimate the error term, we have the sum 

T r /rr\T 


E_ d< = Td - 1 + 


t=i 


(d-2) 2 


by lemma 5.2. Since d ^ 3, we can be sure that 

T—l — T{d — 1) < 0, 
so 

~ / w , d—1 (d — 1) T 

E w ^)(rf-i) f < k ; 


t=i 


(d — 2) 2 T 2 


(5.8) 


(5.9) 


(5.10) 


The first term of (5.7) is easy to bound: since TrOp(/) 2 ^ 2 Bn 2 if |/&| ^ k, we have 


TrOp (/) 2 < k 2 

2 BT " T' 


The final step needed is to bound 

T 


5> T (t)(/,MV)c-, 


(5.11) 


(5.12) 


t=l 
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where, by theorem 4.7, 


|</,M*/) C2S I ^ ||/|| c 2 b||mV|| c2 b ^ K\\f\\l 2B t 


d — l — P 
d- 1 


(5.13) 


for some constant Ah Since WT { t ) ^ T 1 for all t , we can estimate 


fd-l-pV 1 ^ fd-l-p 

hrr hE 1 ' 


t=l 


t=l 


d- 1 


1 (W) 
T (!-^) 2 
1 (d- l)(d- 1 -/3) 
T J 2 : 


(5.14) 


making use of lemma 5.1. As ||/||^ 2 b h 25 k 2 , we end up with 


Y. ™ T ( t )( f , M*/) C 2B = 0 


t=l 


Bk 2 


TP 


2 / ’ 


(5.15) 


and combining with the other two bounds: 


y(/ , B ) = o|^Uo' V( ^ 1)T|Cfl - 2Tl 


T/A 


BT 2 


(5.16) 

□ 


Remark 5.3. By applying the Holder inequality we may derive the alternative bound to 

(5-15): 

j^wrmfMf) a* = °( t1 -i^ 1+1/ p ) > p ’ q > ^ + ^ = 1, (5 - 17) 

which may be useful to prove quantum ergodicity in cases of families of graphs for which 
P —>• 0 in a slow way as B —>• oo. 


6 Examples 

In this section we provide two examples of families of graphs for which our results prove 
quantum ergodicity when quantised with equi-transmitting scattering matrices. 

6.1 Random regular graphs 

There exist several models for chosing a regular graph on n vertices at random [31]. We 
shall consider the set Q n d of simple, d-regular graphs on n labelled vertices. It follows 
from [37, 38] that the size of G n .d obeys 

/ ud d \ 2 

\G n ,d\ ~ J > as oo. (6.1) 
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We make Q n ^d into an ensemble of random graphs by assigning uniform probability to 
each element [32, 31]. These graphs are bipartite with probability o(l) as n —» oo, and 
connected (a fortiori d-connected) with probability 1 — o(l) [39]. 

In order to use theorem 2.2 to prove that quantum ergodicity holds with probability 
1 — o(l), we collect together some prior results showing that such random graphs are 
expanders, and that they do not have too many short cycles. 

Random regular graphs are known to be almost Ramanujan, due to a result of Fried¬ 
mann [40], which had been conjectured by Alon [27]: for any e > 0, with probability 
1 — o(l) a random d-regular graph has all non-trivial eigenvalues /q of its connectivity 
matrix bounded by |//.;| ^ 2 y/d — 1 + e. This means that we can take any j3 < d — 2 yjd — 1 
in theorem 2.2. A weaker bound valid for d even, with a simpler proof, has been given in 
[41], that would also serve our purpose for those values of d. 

For the number of short cycles in a random d-regular graph, such questions have been 
considered in [42]. For d fixed, theorem 4 of [42] reads: 

Theorem 6.1. Let k = k(n) ^ 3 satisfy k(d — l) k ~ l = o(n). Let S = S(n) = 20Ak(d— l) fe 
with A = A(n) > c for some constant c > 1. The probability that the random d-regular 
graph on n vertices has exactly S edges which lie on cycles of length at most k is less than 


e -5(d-i)* Q y /4k 


( 6 . 2 ) 


To apply theorem 6.1 to our situation, let k = | log^x n, and So = L42n 3 / 5 log^x ?zj 
and define A(n) by 

So = 12 A(n)n 3 ^ 5 log^. j n. (6.3) 

Then, for n sufficiently large 

3 < A(n) < 3.5, (6.4) 

and 

k(d - 1) fc_1 = = o(n), (6.5) 

so that the conditions of theorem 6.1 are satisfied. 

The probability that a random d-regular graph has at least So edges which lie on cycles 
of length at most k is, according to theorem 6.1, not more than 


e -^ d ~ 1)k Q y /4k <: 


S=S 0 


= e 


—5 n 3 / 5 


( 6 . 6 ) 


for n sufficiently large, since e/A < 1. 

Let 

T = 771 lo grf-i n = 777 l °Sd -1 B - — logrf.! 
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10 


10 


Then (6.6) implies that 


P ( \Cb, 2 t\ < 42 ( — J log d _! 


> 1 -e 


-5(2 B/d) 3 / 5 


(6.7) 


( 6 . 8 ) 
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that is, with extremely high probability. We can then say that 
'( d-l) T 


P 


BT 2 


-\Cb, 


2 T\ 


0 




id-If 

BT 2 


■\Cb,2t\ ^ 


1400 / 2 \ 9/10 


1 


= P(|C B , 2T K42(^ 3/5 log d _ 1 


d) log d _ 1 (2S/d) 

?)) 


B~ 1/1° 


> 1 — e 


_ —5(2 B/d) 3 / 5 


(6.9) 


by (6.8). 

So with T given by (6.7), and with probability 1 — o(l), the right-hand side of (2.16) 
converges to 0 as B —>• oo, leading to quantum ergodicity for a sequence of random d- 
regular graphs quantised with equi-transmitting scattering matrices. 


6.2 Ramanujan graphs with large girth 

As we have stated in section 2.2, the largest theoretical value for which the parameter (3 
can be taken in theorem 2.2 is (3 = d — 2 yfd--T, and such graphs are called Ramanujan. 

Infinte families of Ramanujan graphs with B oo have been constructed for certain 
values of d only: for d = 3 in [43], for d = p + 1 where p is an odd prime in [44, 45] and, 
more generally, for d any prime power in [30]. An existence proof for bipartite Ramanujan 
graphs for all values of d has recently been given in [46]. 

We will take d = q + 1, where q is a prime power. It is known that equi-transmitting 
scattering matrices of size d x d do exist. In [30] a method of construcing non-bipartite, 
connected, d-regular, Ramanujan graphs on n vertices for a growing sequence of ns is 
given. Furthermore, it is proved that the girth g{B ) of such graphs satisfies 

g(B) ^ ^ log d _! n = | log d _! (J^j . (6.10) 

This girth estimate shows that these graphs are close to extremal, since the Moore bound 
[47, Ch. 23] gives a theoretical upper bound of 

21og d _ 1 n(l + o(l)), (6.11) 

for the girth of a d-regular graph on n vertices. In the case that q is a prime, the upper 
bound (6.10) has been shown to be an asymptotic equality [48]. 

Since the girth g{B) —>• oo as B —» oo, we can take any T < ^g(B ) in theorem 2.2, 
and the Ramanujan property shows that the first term on the right-hand side of (2.16) 
can be made arbitrarily small as B —>• oo. With this choice of T, however, it is also 
clear that |C B) 2 t| = 0, so the second term on the right-hand side of (2.16) is absent, 
proving quantum ergodicity for these d-regular Ramanujan graphs quantised with equi- 
transmitting scattering matrices. 

We end with two remarks concerning quantum ergodicity for equi-transmitting Ra¬ 
manujan graphs: if the bond lengths of the quantum graph are linearly independent over 
Q then we can push T up to g{B) — e. Secondly, the effective logarithmic estimate for the 
decay of quantum variance 


1 








is of the same order as can be rigorously proved in other systems [49, 50, 18, 33], but the 
calculations performed in [7] suggest that the true rate of decay should be algebraic (which 
is also consistent with what is conjectured for more general chaotic quantum systems 
[51, 52]). A decay rate of 1/B has been proved for the quantum variance for the graphs 
studied in [5], but this result aside, going rigorously beyond the logarithmic barrier (6.12) 
seems to be for quantum graphs, as with other systems, a difficult problem. 
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A Matrix inequality 

Let T e N. We recall that we defined 


<T, 


u>r(t):=j T )’ 

I 0, otherwise. 


(A.l) 


It is elementary to calculate the inverse Fourier transform of u>T{t), showing that 
r T 


' —T 


w T {t)e 2nitx d t = w T (x) := ( ? t° S x? X ) > (A.2) 

(1, x = 0. 


For our purposes it suffices to note that wt is everywhere non-negative and wt( 0) = 1, 
for all T. 

The following lemma is taken from [33], p. 1463: 

Lemma A.l. Let U be an N x N unitary matrix, {uj^ =1 be an arbitrary orthonormal 
basis ofU, and A be an N x N matrix. Then if T € N, we have 

N T 

L^2 \ {ujjAuj )E wrWTriA^AU 


—t\ 


N 


(A.3) 


j =1 t=—T 

Proof. Let us denote by 9j the eigenphases of U, so that 


Uuj = e 27Tiej Uj. 


(A.4) 


By expanding the trace, we can write 


N 


T~t 


Uj/c.N 


Tr {A*U t AU~ t ) = E (u j ,A*U t AU~ 

3 =1 
N 

= J2e- 2 ^ t (Au J ,U t Au j ) CN . 

3 =1 


(A.5) 
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By inserting the representation 


N 

Au i = ^2( u k,^-Uj)cNUk, (A.6) 

k =1 

we get 

N 

Tr (A*U t AU~ t ) = ^ m ~ ej)t \{u k ,A Uj ) C N\ 2 . (A.7) 

j,k =1 

We multiply (A.7) by wx(t) and sum over all t, invoking the Poisson summation formula 
to get 

T oo N 

Y2 wr(t) Ti{A*U t AU~ t ) = ^ ^ wt(u + — flfc)|(ufc, Auj) C jy| 2 (A.8) 

t=—T n=—oo j,k=l 

N 

^ wRO) E |(u i ,^'u j ) c jv| 2 , (A.9) 

i=i 

retaining the j = k and n = 0 terms of the sums only. □ 
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